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ABSTRACT 

We examine the fine structure lines ^Pi -^^Pq (492 GHz) and ^^Pi 
(809 GHz) of neutral atomic carbon as bulk molecular gas mass tracers and 
find tliat they can be good and on many occasions better than ^^CO tran- 
sitions, especially at high redshifts. The notion of CI emission as an H2 gas 
mass tracer challenges the long-held view of its distribution over only a rel- 
atively narrow layer in the CH/CI/CO transition zone in FUV-illuminated 
molecular clouds. Past observations have indeed consistently pointed towards 
a more extended CI distribution but it was only recently, with the advent 
of large scale imaging of its "^Pi —^'^Pq transition, that its surprising ubiq- 
uity in molecular clouds has been fully revealed. In the present work we show 
that under typical ISM conditions such an ubiquity is inevitable because of 
well known dynamic and non-equilibrium chemistry processes maintaining a 
significant [C]/[^^CO] abundance throughout Giant Molecular Clouds during 
their lifetime. These processes are more intense in star-forming environments 
where a larger ambient cosmic ray flux will also play an important role in 
boosting [C]/[^^CO]. The resulting CI lines can be bright and effective H2 
mass tracers especially for diffuse (~ 10^ — 10'^ cm^"^) gas while in UV-intense 
and/or metal-poor environments their H2-tracing capability diminishes be- 
cause of large scale CII production but nevertheless remains superior to that 
of ^^CO. The best place to take full advantage of CI's capacity to trace H2 
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is not in the low-z Universe, where large atmospheric absorption at 492 and 
809 GHz precludes routine observations, but at high redshifts (z > 1). 

Key words: Galaxies - starburst - galaxies: sub-mm ~ galaxies: Ly-break - 
ISM: molecules - ISM: abundances - ISM: atoms 

1 INTRODUCTION 

The use of ^^CO and its isotopologues rotational lines to measure H2 gas mass and phys- 
ical properties in galaxies is now a well-established technique (e.g. Dickman et al. 1986; 
Elmegreen 1989; Bryant & Scoville 1996) employed successfully throughout the local Uni- 
verse (z< 0.3) (e.g. Sanders, Scoville & Soifer 1991; Young & Scoville 1991; Solomon et al. 
1992; Downes & Solomon 1998) and at increasingly high redshifts (e.g. Brown & Vanden 
Bout 1991; Downes et al. 1992; Omont et al. 1996; Frayer et al. 1998). In high-z observations 
usually only high-J CO lines J + 1 — > J, (J + 1 ^ 4) are accessible for evaluating the mass 
and physical conditions of H2 gas, thus omitting cooler (Tkm < 30 K) and/or sub-thermally 
excited gas (n < 10'^ cm~^). Recently observations of the lowest ^^CO J=l-0, 2-1 transitions 
in several high redshift objects have been made possible with the Very Large Array (e.g. 
Carilli, Menten, & Yun 1999; Carilli et al. 2002a, 2002b) and in some cases thermalized 
transitions up to fairly high rotational levels (e.g. J=5-4) have been observed. However, 
differential lensing can complicate standard interpretation of line ratios by preferentially 
amplifying the emission of the usually more spatially confined warm and dense gas (e.g. 
Blain 1999) yielding the false impression that such gas properties (and the thermalization 
of high-J CO transitions) characterize the bulk of the molecular gas in such systems. In one 
case of an unlensed high-z starburst large scale sub-thermal excitation of ^^CO J=5-4 has 
been inferred (Papadopoulos & Ivison 2002). 

Another limitation of the standard method is encountered in tracing metal-poor H2 gas, 
which becomes more severe in enhanced FUV radiation fields: the reduced dust-shielding and 
CO self-shielding allows UV photons to dissociate CO while leaving the mostly self-shielding 
H2 intact (Maloney & Black 1988; Israel 1988, 1997). This is now corroborated by observa- 
tions of the metal-poor outer parts of typical spirals (Nakai & Kuno 1995; Arimoto, Sofue, 
& Tsujimoto 1996), globally metal-poor objects like the Magellanic Clouds and Magellanic 
Irregulars (e.g. Madden et al. 1997), and blue compact dwarf galaxies (e.g. Barone et al. 
2000). Unfortunately both types of the aforementioned bias can hinder efforts to measure H2 
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gas mass at high redshifts since high-J ^^CO hnes are those used most often, and significant 
amounts of metal-poor gas reside in UV-intense environments (e.g. in Ly-break galaxies). 

The CI J=l-0 line was one of the first Photodissociation Region (PDR) lines to be 
detected from molecular clouds (Phillips & Huggins 1981), yet its H2 tracing potential has 
gone unnoticed. There are three main reasons for that namely, a) the low atmospheric 
transmission at its rest-frame frequency of ~ 492 GHz precluded, until recently, routine 
observations of molecular clouds in this transition, b) instrumentation constraints in terms of 
receiver sensitivity and the ability to map large areas fast (no multi-beam receivers available), 
and c) an early theoretical prejudice of one- dimensional steady-state PDR models that 
pictured CI distributed only in a relatively narrow CII/CI/CO transition zone on the surface 
of FUV-illuminated molecular clouds (e.g. Tielens & Hollenbach 1985a,b) and thus unable 
to trace the bulk of their mass. Interestingly this prejudice was at times reinforced by 
the instrumentation and atmospheric restrictions which allowed only limited mapping of 
CI emission, which in turn was mainly conducted in regions where theory predicted CI 
should be. These early views on the CI distribution in PDRs could explain some general 
features of its emission such as the surprising robustness of the CI J=l-0 line brightness 
over a wide range of FUV intensities and metallicities (e.g. Kaufman 1999). 

This simple picture is now challenged by two large scale CI surveys of the Orion A 
and B molecular clouds (Ikeda et al. 2002) and the Galactic Center (Ojha et al. 2001) 
that find CI, ^^CO, and ^^CO J=l-0 emission to be fully co-extensive, their intensities 
tightly correlated, and with a surprisingly constant N(CI)/N(^^CO) ~ 0.1 — 0.2 ratio over 
dramatic changes in ambient FUV field and physical conditions. These observations along 
with a considerable body of similar past observational evidence (see Keene 1997 for a review) 
point to a CI distribution fully concomitant with that of CO. The present work uses time- 
varying chemical models under typical conditions found in Giant Molecular Clouds (GMCs) 
to further elucidate this issue, and then examines the potential of CI line emission as an II2 
mass tracer. We focus mainly on the following issues, 

1) The current state of the art of CI observations revealing its distribution throughout molec- 
ular clouds (thus allowing it to be a tracer of their mass), and the role of non-equilibrium 
and/or dynamic processes as well as cosmic rays under typical ISM conditions. 

2) The expected CI line flux densities over a range of physical conditions and their relative 
observational advantage (if any) with respect to the ^^CO transitions. 

3) The uncertainties of H2 mass estimates from CI emission with respect to those associated 
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with the standard method employing CO hnes, as well as the observational prospects of 
using the CI lines to trace molecular gas at high redshifts. 
Throughout this work, we adopt a qo = 0.5 cosmology with Qjy = 0. 

2 THE DISTRIBUTION OF CI IN GMCS 

In an early attempt to account for the mounting number of observations finding CI to be 
more extended than predicted by the models (e.g. Keene et al. 1985) low density PDRs 
(n ~ 100 cm"^) have been advocated. There the CI abundance in the CII/CI/CO transition 
zone becomes larger ([C]/[^^CO] ~ 1 — 100) than in denser ones (n ^ 10^ cm~^) because a 
different, Cl-rich, chemical route is favored, allowing CI to dominate over a larger region of 
the cloud, Av ~ 3 — 7 (Hollenbach, Takahashi, & Tielens 1991). More realistic inhomogeneous 
PDR models showed a deeper FUV penetration and a more pervasive CI emission "coating" 
the clumps making up the GMCs (Meixner & Tielens 1993; Spaans & van Dishoeck 1997). 

The aforementioned schemes can produce a spatial correlation between CI emission and 
that of H2 mass tracers like the ^^CO, ^^CO J=2-l, 1-0 transitions but cannot readily explain 
their tight intensity correlation (Keene et al. 1997). The latter is surprising because e.g. a 
CI/^'^CO intensity ratio should depend on quantities like the surface-layer to total cloud 
mass ratio, the different ambient excitation conditions of the two species (even in clumped 
PDRs CI and ^^CO still reside in different volumes), all expected to vary strongly within a 
molecular cloud and between clouds. The discovery of significant CI emission deep in UV- 
protected sites of cold dark clouds in the Galaxy (e.g. Oka et al. 2001), and M 31 (Israel, 
Tilanus, & Baas 1998) makes the standard steady-state PDR interpetation (with typical 
cosmic ray fluxes) even more difficult. In the few cases where CI J=2-l is also detected, the 
CI (2-l)/(l-0) ratio yields (in LTE) similar gas temperatures to those deduced solely from 
the ^^CO lines or FIR/sub-mm dust continuum in either starburst or quiescent environments 
(Stutzki et al. 1997 and references therein). Finally the striking similarity of CI and ^^CO line 
profiles in the latest large scale survey by Ikeda et al. (2002) suggests a spatial and intensity 
correlation continuing to still smaller scales than those probed by the nominal resolution 
of the particular observations of ~ 0.3 pc (assuming macroscopic motions of much smaller, 
radiatively-decoupled "cells" responsible for line-formation in GMCs, see e.g. Tauber 1996). 
The simplest explanation of all the aforementioned observational results is that CI and CO 
are fully concomitant and trace the same H2 gas mass. 
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2.1 Is the [C]/[CO] equilibrium reachable? 

Non-equilibrium chemical processes can maintain a high [C1]/[^^C0] in cloud regions that 
have not attained chemical equilibrium (e.g. de Boisanger & Chieze 1991; Lee et al. 1996; 
Storzer, Stutzki, & Sternberg 1997), or have done so under turbulent diffusive "mixing" of 
a Cl-rich cloud envelope with its deeper regions (Xie, Allen, & Langer 1995; Xie 1997). The 
most recent investigation of the time dependance of the [C]/[^^CO] ratio and the resulting CI 
line intensities by Storzer et al. examines the effects of abrupt changes of a strong local FUV 
field intensity (Go ~ 10^) caused by cloud-cloud shadowing over timescales of < 10^ yrs. 

In the present work we assume ambient conditions expected for the bulk of the molecular 
gas in GMCs rather than in their few O, B star-illuminated hot/dense spots. These conditions 
consists of: a) initially atomic clouds (the likely precursors of GMCs), b) softer ambient FUV 
fields, and c) a larger range of gas densities. Moreover, processes other than cloud-cloud 
shadowing govern the evolution of GMCs and their ambient conditions, with timescales 
short enough to allow them to maintain a high [C]/[^^CO] ratio throughout the volume of a 
typical CMC during its lifetime. 

In the case of chemical processes the II2 formation timescale provides a good approxi- 
mation for an overall equilibrium to be reached in a molecular cloud that is initially atomic 
(e.g. HoUenbach & Tielens 1999). This timescale is 



where (n) is the average H density and Rf = 3 x 10^^^(T/K)^/^ cm^ s~^ is the canonical H2 
formation rate on dust grains (e.g. Jura 1975). For (n) ~ (20—60) cm~^ and T ~ (70 — 120) K 
(typical for the Cold Neutral Medium HI out of which GMCs form), it is tch ~ (1"~3) x 10^ yrs. 
On the other hand, the dynamic evolution timescale is 



(spherical cloud assumed). This is comparable to tch and is in accordance with the well known 
fact that timescales of ~ (10^ — 10'') yrs characterize a wide variety of processes that fully 
disrupt or otherwise drastically alter typical GMCs. Some of the most important ones are star 
formation with the disruptive effects of O, B, star clusters (e.g. Bash, Green, & Peters 1977), 
turbulent dissipation (Stone, Ostriker, & Gammie 1998; MacLow et al. 1998), and inter- 
cloud clump-clump collisions (Blitz & Shu 1980). The fact that these timescales are shorter 
or comparable to tch suggests that molecular clouds may never reach chemical equilibrium. 
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Driven by the evolution of continuously forming stellar populations throughout a galaxy, 
the average FUV field regulating the CII/CI/CO transition zone in stationary PDR models 
also varies. Recent work by Parravano, HoUenbach & McKee (2003) has shown 



to be the typical time needed for a significant perturbation of the FUV energy density to 



the star formation rate normalized to its Solar neighborhood value. For e.g. U = 10 x Uo it 
is tpuv ~ 5 X 10 yrs, comparable to tch even in quiescent environments (Asfr = !)• 

In order to track the evolution of the [C]/[^^CO] ratio during the aforementioned timescales 
we present results from our homogeneous, time-dependent PDR models in Figure 1. The 
chemistry of gas is represented over a semi-infinite slab, taking into account all the heating 
and cooling processes (see Papadopoulos, Thi, & Viti 2002 and references therein). We have 
also included the latest advancements in coUisional/reaction rates and chemical networks 
which now include most charge-exchange reactions. The latter are particularly important 
in the proper estimate of the [C]/[^^CO] values in heavily shielded, cosmic ray-dominated 
environments with various metallicities (sections 2.2, 2.3). 

The diagrams in Figure 1 confirm that the [CI]/[^^CO] equilibrium timescale and tdyn 
(Equation 2) remain comparable throughout the density range examined. Hence, since self- 
gravity is the plausible main "driver" behind many processes altering the physical conditions 
of a GMC (e.g. star formation; Elmegreen 2002), it can be argued that [CI]/[^^CO] remains 
significant, never settling to its (small) equilibrium values throughout its lifetime. In more 
active star forming environments all processes that disrupt GMCs or alter their ambient con- 
ditions operate on still shorter timescales. The effect of higher star formation rates is easily 
highlighted by Equation 3 since for AgpR > 1? tpuv simply becomes proportionally shorter. 

In the absence of a generally accepted GMC formation and destruction theory the most 
general way to view the effects of elevated star formation on chemical equilibrium is by 
assuming that star formation powers the cycling of interstellar gas between the Warm, 
Cold HI (WNM, CNM), and the molecular gas phases. Then the gas condensation rate 
co'(nH) = d(lognH)/dt = 0sFR/g(nH), where g(nH) = 9M/9(lognH) is the mass fraction per 
logarithmic density interval (Lioure & Chieze 1990). Clearly, irrespective of the particular 
physical processes defining co'(nH), a higher 0sfr induces a faster evolution of any given gas 




(3) 



occur, where k = 1 — 



2, U is the FUV energy density, Uo = 10 -"^"^ergs cm ^, and Asfr is 
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n = 2 10^ cm'^ 




Av 

n = 2 1 0^ cm'^ 




Figure 1. Results from a time-dependent PDR code (see text) for an initially atomic cloud with one-sided illumination by 
a unidirectional FUV field and densities shown at the top of each frame. The assumed metallicity is solar (Z=l) and the 
incident radiation field is Go = 1/2, (corresponding to one-sided illumination by a Habing field). The well-known [C]/[^^CO] 
enhancement in low density gas with respect to denser gas is clearly discernible 

"parcel" across the entire mass-density spectrum of the neutral ISM, leaving less time for it 
to reach chemical equilibrium at any one phase. 



2.2 The role of turbulence 

A relatively recent observational result that is particularly hard to reconcile with steady- 
state models was obtained by the extensive CI J=l-0 and ^^CO J=4-3 imaging of the Carina 
molecular clouds where a close spatial association of their emission throughout the cloud was 
found, suggesting that the same gas emits both lines (Zhang et al. 2001). Low density PDRs 
(~ 10^ cm~^) would be Cl-rich and bright in CI J=l-0 but not in the ^^CO J=4-3 transition 
that has n43 ~ 2 x 10^ cm^'^, while in the bi-stability scenario (section 2.2) a locally much 
higher CR flux (so that n(CR) > 10'^ cm~^) would be needed to explain such association. 

Turbulent diffusion can smooth any initial [CI]/[^^CO] abundance gradient so that it 
becomes almost constant throughout a typical cloud, and in the Carinae molecular complex 
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the clouds are found to be warmer than typical ones in the disk, suggestive of turbulent 
velocity field dissipation (Zhang et al. 2001). The effects of turbulent diffusion on elemental 
abundances have been outlined by Xie, Allen, & Langer (1995) (see also Xie 1997 for a 
review) , and here we will only point out the effects in the kinematically violent environments 
found in e.g. galactic centers and merger systems. In such regions molecular clouds have 
very high velocity dispersions, typically ay ~ 50 km s~^ for the Galactic Center (e.g. Oka 
et al. 1996), that can reach up to ~ (50 — 100) km s~^ in the extreme ULIRG merger 
systems (albeit model-dependent, e.g. Downes & Solomon 1998). These are much higher than 
~ (1 — 5) km s"^ found for GMCs in the kinematically quiescent spiral disks. These large line 
widths in kinematically violent environment correspond to a turbulent diffusion coefficient 
K = (VtL) ~ (10^'^ — 10^^) cm^ s~^, (assuming ay ~ Vt and an invariant characteristic length 
scale L). This is 1-2 orders of magnitude larger than in kinematically quiescent environments 
and can easily yield a constant [C]/[^^CO]~ 0.1 — 1 ratio throughout the cloud (Xie et al. 
1995). The diffusion timescales are short, e.g. for quiescent GMCs tdifr ~ 3 x 10^ yrs, which 
is shortened to ~ 3 x (10"^ — 10^) yrs in dynamically active environments (tdifr oc 1/K). 

Such effects, and a potentially large cosmic ray flux (section 2.3), can easily account for 
the significant [CI]/[^^CO] ratios found in dark and dense clouds, shielded from the ambient 
FUV radiation by very large extinctions (e.g. Ikeda et al. 2002), as well as the similarity of 
CI J=l-0 and HCO+ J=l-0, 3-2, CS J=2-l line profiles (Stark et al. 1997; Johansson et al. 
1994). The latter are particularly puzzling since in shielded environments HCO+ is a product 
of dense, dark cloud chemistry, and its J=3-2 transition has n^ ~ 10^ cm~^ which, along 
with CS J=2-l (ucr ~ 10^ cm~^), imply very dense gas for which standard PDR models 
predict ^^CO to be the dominant carbon-bearing species ([CI]/[^^CO] < 10^'^). Turbulence 
has the strength and the time to redistribute any Cl-rich, PDR-origin gas phase initially 
located in a cloud envelope throughout its denser and UV-shielded interior within its lifetime. 

2.3 Processes controlling [CI]/[^^CO]: the role of cosmic rays 

The previous brief exposition of time-dependent chemical and dynamic phenomena suggests 
that except in few of their densest, most UV-protected and kinematically quiescent regions, 
GMCs may never attain chemical equilibrium while turbulence is capable of "mixing" the 
resulting Cl-rich gas phase throughout their volume, a situation accentuated in star forming 
and/or dynamically active environments. 
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In star forming environments the expected larger cosmic ray flux will also raise the 
[CI]/[^^CO] ratio deep inside a cloud's volume, and several observations have indeed found 
a systematically larger [CI]/[^^CO] ratio in starburst nuclei (e.g. Harrison et al. 1995; Israel 
& Baas 2002). In the starburst M 82, bright CI J=l-0 emission with N(CI)/N(i2C0) ~ 0.5 
is measured over the bulk of its H2 gas (Schilke et al. 1993; White et al. 1994) while the CR 
flux in its central regions is Fcr ~ (170 — 500) x F^^'^ (Suchkov, Allen, & Heckman 1993). 

In Figure 2 the effect of Fcr on the CI distribution is presented, and the CR-induced 
[CI]/[^^CO] enhancement throughout a cloud becomes apparent (see also Flower et al. 1994 
for early such work). Indeed in the dark cloud conditions prevailing deep into the inner 
regions of GMCs a pure PDR origin of CI becomes unattainable but cosmic rays, unlike 
ultraviolet photons, penetrate deep into clouds and at high Ay control both thermal balance 
and chemistry. Such a chemistry can still yield a steady-state Cl-rich gas phase, provided 
n < n(cR). This "turnover" density increases for higher Fcr and marks a CR-controlled 
bi-stability of the chemical network set by the abrupt transition between two possible ISM 
ionization states (des Forets, Roueff, & Flower 1992; Le Bourlot, des Forets, & Roueff 1993; 
Flower et al. 1994). For n < n(cR) the gas has a high ionization fraction and thus a high 
dissociative recombination, allowing 11+ to be the dominant ion and initiate a set of 
chemical reactions that yield [CI]/[^^CO] ~ 0.1 — 0.2. For n > n(cR), Hg" becomes dominant 
and a different chemical network yields a much lower [C]/[^^CO] < 10"'^. 

For the Galaxy it is n(cR) ~ 5 x 10'^ cm~^ (Schilke et al. 1993 and references therein), 
encompassing a large fraction of the molecular gas in GMCs. For totally shielded environ- 
ments (Go = 0) and Fcr = f'§^^^ our models yield [C]/[^^CO] > 3 for n ^ 10^ cm-^, and 
[C]/[^^CO] ^ 0.02 for n ^ 5 X 10'^ cm~'^. A ten-fold increase in Fcr boosts the aforemen- 
tioned densities by an order of magnitude, for similar [C]/[^^CO] values. These values are in 
good accordance with those reported in the literature even though controversies remain and 
a major benchmarking effort of the various PDR codes is under way. In an important recent 
development the canonical value of Fcr in quiescent diffuse clouds in the Solar neighborhood 
may have to be revised upwards by a factor of ~ 40 (McCall et al. 2003) thus opening up 
the possibility of a mainly CR-origin of the CI emission observed to be widely distributed 
in the quiescent molecular clouds in the disk of the Galaxy. 

In Figure 3 we explore the same effects but now scale both the CR flux and the FUV 
ambient radiation field, since their average values are expected to rise in tandem in realistic 
star forming environments. In this case enhanced CII production at large FUV intensities 
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n = 2 10''cm^-CR=1 n = 2 10" cm^ - CR=10 n = 2 10" cm ^ - CR=100 




0.1 1.0 0.1 1.0 0.1 1.0 

Av Av Av 



Figure 2. Abundance distributions at t = lO'^ yrs, for identical initial conditions as in Figure 1. The values of Fcr are indicated 
in the panels, normalized to the Galactic value. The ambient Go = 1/2, corresponding to one-sided cloud illumination by a 
Habing field. In most cases [CI]/[CII]> 1 for Av ^ 1. The density range chosen characterizes the bulk of H2 in typical GMCs. 

limits the cloud volume where CI is dominant, but except for diffuse gas irradiated by strong 
fields, CI remains the dominant form of carbon for Av > 2. 

2.4 CI distribution in metal-poor environments 

The expected evolution of metallicity in galaxies with an overall trend of lower values at 
earlier epochs is now verified by numerous observations (e.g. Mehlert et al. 2002). A trend 
of higher star formation rates in galaxies at higher redshifts is also deduced, although the 
details are far from clear (e.g. Ellis 2001; Somerville, Primack, & Faber 2001). It is therefore 
anticipated that galaxies with low metallicity and elevated star formation rates (and thus 
FUV and CR fiux), will be common at high redshifts. Ly-break galaxies found at z ~ 3 
with moderate star formation rates of ~ (5 — 10) MQyr^^ (Steidel et al. 1996), and sub-solar 
metallicities (e.g. Pettini et al. 2000) are typical such objects. 

In low-met allicity environments with pervasive FUV radiation powered by star forma- 
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n = 2 10" cm ^ - CR=1 - Go=0.5 n = 2 10" cm ^ - CR=10 - Go=5 n = 2 10" cm ^ - CR=100 - Go=50 




Figure 3. Abundance distributions at t = lO'^ yrs, assuming the same initial conditions as in Figure 1. The values of Fqyi a^re 
indicated in the panels, normalized to the Galactic value. The ambient Go is now scaled proportionally but also multiplied by 
1/2 corresponding to one-sided cloud illumination. 



tion, the enhanced CO photodissociation diminishes its spatial extent per molecular cloud 
and thus its H2-tracing capability (Pak et al. 1998; Bolatto, Jackson, & Ingalls 1999). Ob- 
servational work using the AST/RO facility in the South Pole Station confirms these expec- 
tations by finding enhanced CI/^^CO J=l-0 line ratios in the LMC and SMC (Stark et al. 
1997; Bolatto et al. 2000a). These data, albeit sparse and with a large dispersion, suggest 
Ici/Ico oc Z~^/^, supported also by observations in the low-metallicity star forming irregular 
IC 10 (Bolatto et al. 2000b and references therein). 

In Figure 4 the effects of such conditions on the CI distribution in solar and sub-solar 
metallicity environments are contrasted, and it can be seen that CI remains the dominant 
form of carbon for the bulk of the cloud (Av ^ 2) even in the metal-poor case. Thus in the 
metal-poor and modestly FUV- enhanced environments, characterizing many galaxies at high 
redshifts, CI remains a good tracer of molecular gas mass. 

A large change is seen however in the distribution of CII that now, unlike in the metal- 
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z=0.25 n=2 1 0' cm ' CR=1 Go=5 z=1 .0 
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X(CI) 

X(CII) 




10.0 



Figure 4. Carbon species abundance distributions at t = 10^ yrs relative to the total carbon content. The cloud is assumed 
initially atomic, and Fcr, Go correspond to an underlying star formation rate of ~ IOM0 yr~^, typical of Ly-break galaxies. 
The sub-solar metallicity of Z=0.25 is also typical for Ly-break galaxies (Pettini et al. 2000). 

rich case, it becomes more abundant than CO even deep into the cloud (yet still less than 
CI). Apart from the enhanced FUV photoinization of CI at low Ay's, caused by the reduced 
dust shielding, low metallicity further enhances CII at the deeper parts of such clouds at 
the expense of CO because the main route of CO destruction at high Ay is via the reaction 
path: CO + IIe+ —>■ C^ + O + He. Once formed, ionized carbon will capture an electron 
from a metal atom in the gas phase but at low metallicities the charge exchange reactions 
involving the metallic atoms and ions are less efficient in converting C"*" into C while He^ is 
formed in large abundance by the larger cosmic rays flux and thus a significant part of the 
carbon will be in the form of CII even at large Ay's. 
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3 THE INTENSITY OF THE CI LINES 

The critical densities of the fine structure CI fines ^Pi ^^Pq (z^io = 492.160 GHz) and 
^P2 — s>^Pi (z/21 = 809.343 GHz) are uio ~ 500 cm"^ and n2i ~ 10^ cm~^, comparable to 
tfiose of tfie widely observed ^^CO J=l-0, 2-1 transitions, wfiile tfieir upper energy levels 
are Ei/k = 23.6 K, and E2/k = 62.4 K. Tfie velocity-integrated flux density of a line is 



wfiere Dl = 2cH^^(l + z — + z) is tfie luminosity distance for a redsfiift z, Lx is tfie 
area/velocity-integrated line brigfitness temperature (in K km sec~^ pc^), and its rest- 
frame frequency. Hence, for any of tfie CI J=l-0, 2-1 and tfie ^^CO J+1— i>J lines tfie velocity- 
integrated flux density ratio is 



/^.ScidV /^ei^ RCI ... 

/^^ScodV WJ Rj+i,j' ^ 
wfiere Rci = (Tb(CI))/(Tb(COi_o)), and Rj+i,j = (Tb(J + 1, J))/(Tb(l, 0)) for ^^CO. Tfie 



brigfitness temperature averages are over velocity and emitting area. 

An unresolved (or barely resolved) object, wfiose various ISM lines are accessible by a 
multiplicity of receivers, will fiave an observed CI flux density larger by tfie (z/ci/t'co)^ factor 
relative to tfiat of lower frequency CO transitions. Interferometric scaled-array observations 
'syntfiesizing' tfie same beamsize at widely different frequencies wfiile utilizing tfie same col- 
lecting area maintain tfiis advantage in tfie case of resolved objects. A flux density advantage 
of a CI line may remain even wfien Va/^co ^ 1? simply because CI lines are mucfi easier 
to excite and tfius can maintain tfieir intensities intact for diffuse and cooler gas wfiile tfie 
Rj+i J ratios diminisfi (see also Kaufman 1999). For example tfie ^^CO J=4-3 and J=7-6 
lines fiave critical densities n43 = 2 x 10^ cm~'^, nyg = 3 x 10^ cm~'^ and E4/k ~ 55 K, 
E7/k ~ 155 K, all mucfi fiigfier tfian tfie corresponding quantities of tfie CI J=l-0, 2-1 lines 
emitted at similar frequencies. 

Tfius for any detected ^^CO J+l—>-J (J+1^3) line, bright CI emission may also he 
present, particularly for sub-thermally excited, and/or cooler gas (low Rj+i^j. Tfiis is of 
particular importance for tfie mounting number of fiigfi redsfiift (z>2) CO detections, wfiere 
only tfie J + 1 ^ 3 lines are usually observed wfiile for sucfi redsfiifts tfie CI lines no longer 
undergo large atmospfieric absorption. 

In Table 1 we summarize tfie typical brigfitness temperature line ratios measured in 
various ISM environments in tfie local Universe gleaned from an extensive searcfi in tfie 
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literature. Large atmospheric absorption has prevented any large surveys of ^^CO J+l^J, 
J+l>3 or ^^CO J=3-2, so we use a single-phase Large Velocity Gradient (LVG) code, 
constrained by the CO line ratios available in each case, to estimate expected flux densities 
for J+l>3. The results are then used to calculate detection advantage ratios Da, defined as 
the signal-to-noise (S/N) ratio normalized by that of ^^CO J=l-0 for equal integration time 
and velocity resolution. A general expression for the system sensitivity is 



1/2 

1/2 ' 



sys 



Z/TintAV 



ASsys(z/), (6) 



where = (1 + z)^^!/^ is the observed frequency of the line. Tint is the integration time, 
AV is the velocity resolution, and ASsys(z^) describes the overall system sensitivity and 
incorporates instrumental (e.g. collecting area, receiver temperatures) as well as atmosphere- 
dependent factors (e.g. absorption, signal decorrelation factors in case of interferometers). 
From Equations 4 and 6 the S/N can be expressed as 



5/2 



[T-intAV) 

ASsys(z^) 



1/2 



5 



(7) 



where (Tb) is the velocity /area- aver aged line brightness temperature. For frequencies of ~ 
30 — 250 GHz and future modern instruments like ALMA, located in dry sites, ASsys(z^)/z^^^^ 
changes < 50%^. Thus setting ASsys(z^) oc z/^/^, and for a common Tint and AV, 

Da = = {^J^ioYR (8) 

expresses the detection advantage of any spectral line with rest-frame frequency u^, nor- 
malized to ^^CO J=l-0, with R being their brightness temperature ratio. This assumes 
a bandwidth fully covering both lines under comparison, which will be true for planned 
mm/sub-mm telescopes. 

The Da values for the ^^CO and CI lines are shown in Figure 5 where it is obvious that 
^^CO J=4-3 marks a broad excitation "turnover" beyond which Da(^^CO) starts declining. 
This is not surprising since for most of the gas in GMCs T^in < 50 K, and n(H2) < 10"^ cm~^, 
rendering the ^^CO J+1— s>J, 3 lines sub-thermally excited and underluminous. The wide 
range of Da(CI) values reflects the considerable range of Rci ratios observed in various 
environments, which is not simply correlated to that of the ^^CO line ratios, e.g. a large 
Rj+i^j not necessarily mean also a large Rcj. The CI J=l-0 line has a clear advantage over 
^^CO transitions in quiescent environments and maintains a sizable one in starbursts. The 

1 See |http:// www . alma. nrao . edu /info / sensitivities | Table 1 
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Very Quiescent 



1-0 2-1 3-2 4-3 [CI] 5-4 6-5 7-6 [CI] 8-7 9-8 



Starburst 



■■■1. I 



1-0 2-1 3-2 4-3 [CI] 5-4 6-5 7-6 [CI] 8-7 9-8 



Quiescent 



■■ ■ . 



1-0 2-1 3-2 4-3 [CI] 5-4 6-5 7-6 [CI] 8-7 9-8 



Extreme Starburst 



-III. I 



1-0 2-1 3-2 4-3 [CI] 5-4 6-5 7-6 [CI] 8-7 9-8 



Dynamically Active 



■ll. 



1-0 2-1 3-2 4-3 [CI] 5-4 6-5 7-6 [CI] 8-7 9-8 



Metal-poor starburst 



I 



1-0 2-1 3-2 4-3 [CI] 5-4 6-5 7-6 [CI] 8-7 9-8 



Figure 5. The values of Da for the ISM environments listed in Table 1 (see text). A dispersion of ~ 15% was assumed for all 
ratios (typical error of the mean value reported by surveys). For dynamically active and extreme starburst environments we 
adopt the much larger range of the CO line ratios found around the mean values reported in Table 1. For metal-poor staburst 
we set Z=0.25, and adopt the scaling Rci(Z) = Z~^''^Rci(l) (Bolatto et al. 2000b), all other ratios are assumed the same as 
in starbursts. 



CI J = 2 — 1 is the only line in the molecular ISM that maintains high Da at such high 
frequencies under most conditions, a conclusion subject to the uncertainty of the assumed 
CI(2 — 1)/(1 — 0) ratio which has not been widely measured. For metal-poor starbursts 
Figure 5 can be somewhat misleading since in such environments CII may dominate over 
the bulk of the molecular gas mass with both the ^^CO and CI lines being rather weak. 
Nevertheless in the case of a detected ^^CO line in such an environment, CI line emission is 
expected to be bright as well. 

It must be emphasized that values translate to a practical observational advantage 
provided the relevant lines do not fall near any of the (many) atmospheric absorption win- 
dows within the 30 — 250 GHz frequency range, and that the set of ^^CO, CI lines under com- 
parison is accessible within this frequency interval. The latter implies redshifts of zi_o ^ 1 
and Z2_i > 2.25 so that the CI J=l-0, 2-1 lines fall within the aforementioned range. 
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A single receiver observing various molecular ISM lines redsliifted within its tuning range 
from several objects distributed over a redsliift range is a case more relevant to current 
instrumentation capabilities (limited set of receivers per telescope) and blank field spectral 
line surveys to be conducted with future instruments per observing band (e.g. Blain et al. 
2000). Then Equation 8 is no longer relevant and one has to simply look at Equation 7 and 
the expression in the brackets which changes < 3% for z ~ 1 — 10. This allows ^^CO, CI lines 
to trace similar amounts of gas mass over a progressively larger redshift range as long as 
sufficiently excited lines with progressively larger rest-frame frequencies become accessible. 
This is simply the spectral-line equivalent of the well-known favorable K-correction that 
keeps the FIR/sub-mm dust continuum of galaxies undiminished (and thus detectable) over 
a large range of redshifts. For ^^CO lines this is beautifully demonstrated by ^^CO J=l — 
observations of nearby QSOs that measure H2 masses similar to those found in high redshift 
QSOs and radio galaxies detected through higher-J ^^CO lines (Evans et al. 2001). 

Studies examining the luminosity of molecular and atomic ISM lines in the context 
of emission-line searches at high redshifts (e.g. Blain et al. 2000) do not explore the H2- 
tracing potential of the CI lines. Moreover, by making favorable assumptions about the gas 
excitation they present optimistic predictions of flux densities, particularly for high-J ^^CO 
transitions. This is why some of their brightness temperature line ratios are unrealistic, e.g. 
CI (2 - 1)/(1 - 0) ~ 2.25 (Table 1 in Blain et al. 2000), while the highest value observed to 
date is ~ 1, in M 82 (Stutzki et al. 1997) and the Orion KL region (Yamamoto 2001). 

There is no reason to believe that the variety of environments indicated in Figure 5 will 
not be also present at high redshifts. Even if starburst environments are more frequently 
encountered there, the possibility of large scale sub-thermal excitation as suggested by re- 
cent measurements in the extremely red object (ERO) and starburst HR 10 at z ~ 1.45 
(Greve et al. 2003), cannot be discounted. The excitation turnover at the frequencies of the 
CI J=l-0 and ^^CO J=4-3 lines implies that, with the exception of the CI J=2-l line, the 
gas mass tracing capability of higher frequency lines {u > 492 GHz) is suppressed except in 
star-forming environments. For essentially the same reasons molecular gas mass estimates 
based solely on high-J ^^CO lines observed at high redshifts carry large uncertainties (fac- 
tors of ^10). Summarizing, while metal-rich gas fueling star-formation is bound to figure 
prominently in high-J ^^CO line emission at high redshifts, it will not be the only gas phase 
present and large amounts of colder, diffuse molecular gas may be lingering e.g. "in-between" 
successive star forming events beyond the current star-forming sites. 
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Here we must mention that the case for observing the low ^^CO J=l-0, 2-1 hnes remains 
strong even if the CI hnes can dehver more flux per H2 mass for e.g. cool and diffuse gas 
reservoirs. This is because of the current commissioning of large single dish radio telescopes 
like the Green Bank Telescope (GBT) (tunnable up to 115 GHz), and future instruments as 
the expanded VLA (EVLA), (continuously tunnable below 45 GHz), and a Square Kilometer 
Array (SKA) equipped with high frequency receivers (~ 20 — 45 GHz). These will be able to 
conduct very sensitive observations of the low-J ^^CO transitions at high redshifts (Carilli 
& Blain 2002), bypass the excitation biases mentioned previously, and sensitively trace 
molecular gas in the distant Universe. 



4 DEDUCING H2 MASS FROM CI LINE INTENSITIES: THE 
UNCERTAINTIES 

The CI lines have small to moderate optical depths (rci ~ 0.1 — 1) being, over large scales, 
as optically thin as ^'^CO (Ojha et el. 2001; Ikeda et al. 2002). This property along with their 
simpler 3-level partition function and easy thermalization makes them more straightforward 
H2 gas mass tracers than the usually optically thick ^^CO lines. 
For an optically thin line the emergent brightness is 
hz/ui 



I 



An 



-A„,N 



(9) 



where Nu is the column density of the upper level, and (puem, is the intrinsic line profile. The 
received brightness is I,^ = (1 + z)^H^^^^ (see e.g. Ivison et al. 1996). Thus in terms of the 
usually reported observable I = / S^dV (S,^ = / l,^dfl) the CI beam-averaged column density 
over the telescope beam fib can be expressed as 
_ 4vr(l + z)3lci 

where the Qui = Nu/Nci < 1 excitation factor depends on (n, Tk), and the radiation field 
(see Appendix). In terms of the gas mass contained within f2b the last expression becomes 

where we used Ab = (1 + z)^'^DlQ\^ as the corresponding area of the beam at the redshift of 
the source, Xci = [C]/[H2], and /i = 1.36 corrects for the mass in He. In astrophysical units. 
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where we set Ho = lOOh' km s~^Mpc~^ 

For optically thick emission the latter expression yields a lower limit on the gas mass 
and remains valid for optical depths emerging in radiatively decoupled "cells" much smaller 
than the emitting region, but with Aui replaced by /5uiAui (/5 = (1 — e~'^)/r: photon escape 
probability). In molecular clouds where tqi < 1, an unknown value of Pui introduces an 
uncertainty of < 60%. A greater source of uncertainty stems from the assumed excitation 
conditions that determine Qui. For high redshift objects where only one or two lines are 
usually detected there can be no practical constraints on the gas excitation and the minimal 
assumption of LTE is usually made in order to deduce M(H2). For typical n = (100 — 5 x 
10^) cm-3 and Tk = (20 - 60) K it is Qio/Q£'^^^ ~ 0.35 - 1 and Q2i/Q2i™^ ~ 0.15 - 1. 
These uncertainties can be larger if there are no available temperature constraints from e.g. 
measurements of the dust continuum. Predictably Q21 has the largest variation, making the 
CI J=2-l transition the more uncertain II2 mass estimator of the two CI lines. Nevertheless 
the ^^CO J=7-6 line at a similar frequency (806 GHz) is much more sensitive to ambient 
conditions, and hence an even more uncertain H2 mass tracer. Similarly the CI J=l-0 with 
a smaller Qio variation is tracing H2 gas more reliably than the ^^CO J=4-3 (461 GHz) line 
even though the latter can be significantly stronger. 

In deducing M(H2) from CI lines we assumed CI distributed throughout a typical molec- 
ular cloud but also no significant "contamination" from any other gas phase. This is indeed 
the case since carbon in HI gas is mostly ionized with e.g. [CI]/[CII] ^ 3 x 10^^ in the 
CNM phase (Wolfire et al. 2003). The small CI line optical depths are an advantage when 
estimating H2 mass but they also add Xci to the list of uncertainties. Given the variety 
of processes affecting [CI]/["'^^CO] and thus Xci, it seems that the latter can vary widely. 
However the very same processes also affect Xco = [^^C0]/[H2], which is nevertheless found 
to be relatively constant (~ 10~^) over large scales for clouds with >2 (e.g. van Dishoeck 
1992). This is in accordance with the constant N(CI)/N(12qo) ~ 0.1 - 0.2 found by large 
scale CI, ^^CO, ^^CO J=l-0 surveys (e.g. Ikeda et al 2002) and suggests a dominant physical 
process that along with turbulent diffusion sets a relatively constant [C]/[^^CO] throughout 
a CMC. However, as was the case for ^^CO, establishing average Xci values will probably 
be done observationally, with large surveys conducted over a broad range of conditions and 
discerning any general trends of [C]/[^^CO] (and thus Xci) in active and/or metal-poor 
environments. 
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5 OBSERVATIONAL PROSPECTS AT HIGH REDSHIFTS 

The potential advantage of the CI hnes as H2 mass tracers can be fully realized only at high 
redshifts. For zio > 1 the CI J=l-0 line is redshifted to < 250 GHz where a combination of 
existing large- aperture instruments (single dish and interferometers) and a much improved 
atmospheric transmission enables sensitive observations. For the J=2-l transition the afore- 
mentioned limit is Z21 > 2.25, a typical redshift of the new population of sub-mm bright, 
starburst galaxies (Chapman et al. 2003). The commissioning of interferometers in excellent 
sites, namely SMA on Mauna Kea, Hawaii and ALMA in Llano de Chajnantor, Chile, will 
allow routine sensitive observations at sub-mm wavelengths. In such sites the atmospheric 
transmission even at ~ 345 GHz remains exceptional most of the time (~ 80% at zenith, 
except around the 325 GHz atmospheric absorption feature), and the CI lines become ac- 
cessible to sensitive measurements for zio > 0.4 and Z21 > 1.3. Particularly interesting is the 
case when both CI lines become accessible at lower frequencies, thus allowing an estimate of 
gas temperatures from their intensity ratio (assuming LTE) independently of those based 
on dust continuum or ^^CO transitions. 

Currently all high redshift objects detected in ^^CO are the sites of large scale starbursts 
and are good targets for CI observations, which hold the promise of uncovering H2 gas 
untraced by the high-J ^^CO transitions. Recently an opportunity to drastically increase their 
rather small numbers has emerged with the successful spectroscopic redshift determination 
of several sub-mm bright starbursts at high redshifts (Chapman et al. 2003). Such objects 
are metal-rich and are found to be prominent ^^CO emitters (e.g. Neri et al. 2003). The 
recent detection of weak ^^CO J=3-2 in a Ly-break galaxy at z=2.7 (Baker et al. 2003), 
with its CI J=2-l line redshifted to ~ 218 GHz, offers the possibility of testing the potential 
of CI lines as H2 gas tracers in metal-poor and modestly UV-intense objects. 

The recent detection of massive gaseous disks at high redshifts (Genzel et al. 2003) are 
currently a highly contentious issue since their presence at such early epochs is hard to 
reconcile with standard hierarchical models of galaxy evolution. The much larger spatial 
distribution of the CI line emission with respect to high-J ^^CO transitions in typical disks 
like the Milky Way, holds the promise of yielding a much better picture of disk kinematics 
and total dynamical mass (by encompassing its flat rotation curve part), and thus is an 
excellent alternative. 

There are currently only two reported detections of CI lines at high redshifts, namely in 
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the Cloverleaf quasar (z=2.56) (Barvainis et al. 1997; Weiss et al. 2003), and the IRAS galaxy 
F10214+472 (z=2.28) (Brown & Vanden Bout 1992)^. Both are strongly lensed systems and 
particularities of differential lensing may be responsible for their very different CI/^^CO flux 
density ratios. Clearly CI observations of non-lensed high redshift objects are needed. 

Interestingly in the case of the Cloverleaf, and based on the assumption of a concomitant 
CI emission, Weiss et al. find an II2 gas mass similar to that deduced from the well-sampled 
dust continuum or from the multi-transition ^^CO observations, lending support to the main 
argument of this work. 

Finally, the markedly more difficult CI observations in the local Universe remain valuable 
since large scale CI, ^^CO, and ^'^CO surveys in a variety of environments are needed to 
establish the range of Xci and "cross-calibrate" the two methods of estimating H2 gas mass. 
It is particularly interesting to observe intensely star-forming and/or metal-poor objects in 
order to search for [C]/[^^CO] enhancements and its possible metallicity dependence. 

6 CONCLUSIONS 

In this work we examined the prospects of using the two lines of atomic carbon at 492 GHz 
(J=l-0) and 809 GHz (J=2-l) as alternative H2 mass tracers in galaxies. Our conclusions 
can be summarized as follows: 

1. Non-equilibrium chemical states, turbulent diffusive mixing, and cosmic rays can eas- 
ily make CI ubiquitous throughout a typical Giant Molecular Cloud. This is now firmly 
supported by large scale CI, ^^CO, ^'^CO surveys of such clouds in the Milky Way, which 
also find a constant average [C]/[^^CO] abundance ratio over the bulk of their mass. 

2. In star forming and/or dynamically fast-evolving environments the processes respon- 
sible for the ubiquity of CI in GMCs will further enhance its abundance. In UV-intense and 
metal-poor objects large scale CII production diminishes the H2-tracing capability of CI, 
which nevertheless remains better than that of the now severely dissociated ^^CO. 

3. In terms of emerging fiux density as a function of excitation conditions CI lines have a 
potentially strong advantage with respect to ^^CO lines in tracing H2 gas mass under diffuse 
and/or cooler conditions as well as in the metal poor environments of mild starbursts. 
However the large atmospheric absorption at their rest-frame frequencies allows such an 
advantage to materialize at high redshifts (z > 1). 

^ This detection was under dispute but it has been recently confirmed (Papadopoulos et al. 2004) 
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4. The simpler 3-level structure of the CI hues, their small optical depths and modest 
excitation requirements with respect to ^^CO lines of similar frequency make the associated 
H2 mass estimates more straightforward. However the usually unknown (particularly for high 
redshift sources) mean excitation conditions introduce large uncertainties in such estimates, 
which are reduced considerably if an estimate of the gas temperature is available. 

5. Concentrated observational effort in the local Universe remains critical. Large scale 
CI, ^^CO, and ^^CO surveys of GMCs in a variety of environments will help address issues 
of optical depths, [CI]/[H2] abundances and its variations, and "cross-calibrate" the two H2 
gas mass tracing techniques. 
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Table 1. CO, CI line ratios, typical LVG solutions 



Environment^ 


T3 b 
H2I 








Hci (r2i) 


LVCr solutions 


Keierences 


Very quiescent 


0.4 


0.15 


9 




0.10 (0.22) 


10, 10^, 10-5 (40) 


1, 2, 3, 4 














10, 300, 10-4 (60) 




Quiescent 


0.6 


0.30 


5 




0.15 (0.30) 


20, 300, 10-* (40) 


5, 6, 7, 8 














10, 10^, 3 X 10-5 (60) 




Dynamically active 


0.95 




10 


10 


0.10 (0.22) 


60, 10^, 10-5 (25) 


9, 10, 11, 12 


Starburst 


0.95 


0.65 


15 


15 


0.15-0.30 (0.5) 


30, 300, 3 X 10-" (40)'' 


13, 14, 15, 16, 17 














60, 10^, 10-5 (40) 




Extreme starburst^ 


0.95 


0.65 


25 


25 


0.15-0.30 (0.5) 


70, 10^, 10-5 (60) 


13, 14 



" Very quiescent: M31 cold clouds, Quiescent: Milky Way disk, Dynamically active: Galactic Center, Starburst: various line 

surveys in several objects, Extreme starbursts: ULIRG nuclei 

^ The I^CO (2-1)7(1-0), (3-2)/(l-0) brightness temperature ratios 

The '^■^CO/^^CO brightness temperature ratios for J=l-0, 2-1 
<^ The l^CO/CI J=l-0 and r2i=CI(2-l)/(l-0) brightness temperature ratios 

" LVG solutions: Tk(K), n(H2)(cm-3), Xco/{dV/dr)([km s-ipc-i]-i), ([i^cO/i^CO] = 40,60), Xco = [^^cO/Ha]. In case 
of serious "degeneracies" typical solutions from both parts of the parameter space are displayed (e.g. in starbursts). 
f Best solutions occur for [i^CO/^^CO] = 40 

^ Indicative ratios only, large range of values observed we adopt the R32 and CI values from starbursts. 

1: Israel, Tilanus, & Baas 1998, 2: Allen et al. 3: Fixsen et al. 1999 {COBE, outer Galaxy), 4: Wilson 1997 (M 33 clouds), 5: 
Disk cloud surveys by Chiar et al. 1994; Hasegawa 1997; Sakamoto et al 1997; Falgarone et al. 1998; Sanders et al. 1993 (for 
CO 3-2); Digel, de Geus, & Thaddeus 1994, 6: Disk cloud surveys by Solomon, Scoville, & Sanders 1979; Rickard & Blitz 1985; 
Polk et al. 1988; Falgarone et al. 1998; Wilson et al. 1999, 7: Fixsen et al. 1999 (inner Galaxy); 8: Gerin & Philips 2000; 9: 
Langer & Penzias 1993; 10: Sawada et al. 2001; 11: Ojha et al. 2001; 12: Fixsen et al. 1999 (Galactic center); 13: Braine & 
Combes 1992, Aalto et al. 1995, Papadopoulos &; Seaquist 1998; 14: Devereux et al. 1994; 15: Stutzki et al. 1997; Yamamoto 
et al. 2001; 16: Israel, White, & Baas 1995, Petitpas & Wilson 1998; Israel & Baas 2002, 17: Stark et al. 1997; Bolatto et al. 
2000a, 2000b 
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APPENDIX A: THE 3-LEVEL SYSTEM AND THE CI LINES 

The rate equations for a 3-level system in steady state under the assumption of A20/A21 <^ 1, 
A20/A10 ^ 1, as is the case for the CI fine structure levels, are 



^ = Bi2Ji.2ini - (B21 J,.2i + A21) n2 + [Ci2ni + Co2no - C2in2 - C2on2] = 0, (Al) 



= (B2iJi/2i + A21) n2 + Boi JiyjoUo — (Aio + BioJi^n, + B12J21) ni + 

[C2in2 + Como - Ci2ni - CiqUi] = 0, (A2) 

where J^, = 2hi/^/c^ [exp(hi//kTb) — is the background radiation field (assumed isotropic 
and with no line contribution), Cui = n7ui are the rates of collisional de-excitation (u 1) 
related to those of the reverse process through the principle of detailed balanced Ciu/Cui = 
(gu/gi)exp(— liz/ui/kTk). The term in the square brackets is the net collisional population/de- 
population of the levels, and the rest are the stimulated and spontaneous processes. 

The stimulated excitation and de-excitation coefficients are related as guBui = giBiu 
and to the spontaneous emission coefficient Aui = (2hz/yi/c^)Bui. The latter can be used to 
demonstrate that B20/B21 ^ 1 and B20/B10 ^ 1, also assumed in the previous equations. 
After some reordering and replacing all the Bik factors we get 

(C21 + C20 + f2iA2i) n2 - [C12 + g2/gi (f2i - 1) A21] ni - Co2no = 0, (A3) 

(C21 + f2iA2i) n2 - [Cio + C12 + fioAio + gs/gi (f2i - 1) A21] ni + 

[Coi + gi/go (fio - 1) Aio] no = 0, (A4) 
with f = exp(hi//kTb)/[exp(hi//kTb) — 1]. Equations A3, A4 along with the "closure" relation 
n2 + ni + no = nci, (A5) 
are used to solve for the level populations. We find the following expressions. 



ni/nci ^ „ A, C2o\ nai Coi A , C2o\ A , , n2i' 



(A6) 



K-i ■ - V ' c^j n ■ C02 V C2iy V ^ n 

'C12 + Cio\ Uio 



1 + ^ (fio - 1; 



iO \ ^01 



Coi / n 



© 2004 RAS, MNRAS 000, [T]-?? 



24 Papadopoulos, Thi, and Viti 



n2/nci _ 

1 + - (f21 

and 



n 

'C21 



Clo + C12 g2 

— r + r *^ 2^ 

*^21 gl 



C21 / n C20 



^ /C21 + C20\ 

7^ X 1 + — flO - 1 

V C12 / n go V Coi / 



nio 
n 



(A7) 



no/nci _ Cio A 1121 \ / C20 



l + fio(l + &^l^ 



1 + ^ (f2i - 1 



Cio/ n J Coi 
C20 + C21 ^ n2i 
n 



(A8) 



3I V C12 

It is Qio = ni/iici and Q21 = n2/nci (see main text) while the critical density of a transition 
u — > 1 is n^i = Aui/ (Z]k7uk) (e.g. Jansen 1995), where the sum is over all collisional processes 
out of the level (u). In the literature often only the 7ui coefficient or the sum of just the 
downwards rates is used resulting to an often serious overestimate of Uui- For the 3-level 
atom it is n2i = A2i/(72o + 721) and nio = Aio/(7i2 +7io) (nio > n2i). The expression for K 
is given by 
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In the case n2i/n (and thus also nio/n 0), the last expressions yield 
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The last relations yield 



C 



01 



Coi 
C02 



Cio 
C12 



= 2ii = g^g-E2i/kTk ni ^ Coi ^ gl^_Eio/kTk^ 

ni C21 gl ' no Cio go 

which are the LTE values while the expression for K becomes 



• (A9) 
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(A15) 
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Thus from Equations AlO and A14 it is 

ni 1 g^e"^i/'''^k 

^ ~ 1 + + ~ go + gie-Ei/kT, + g2e-E2/kTk ' 

which its expected LTE value. In the last derivations as well as in all previous computations 
we made use of the permutation relation C02C21C10 = C20C01C12 that results from the 
detailed balance equations. 

In the radiation-dominated limit nio/n 00 (and thus also n2i/n 00) we eventually 
deduce 
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which are the values expected for levels fully regulated by the background radiation field. 

The case where n2i ^ nio allows partial thermalization of the 3-level system since there 
is a density domain where nio/n while n2i/n remains finite and could be even < 1. We 
then find the limit 
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Interestingly this still deviates from its LTE value unless two more conditions are met, 
namely a) f2i ~ 1 (negligible background at the u = V21 frequency), which yields 



ni 
no 



^Ig-Eio/kTk. 



1 + 


(1+ 


C20^ 
C21) 


C02 


1 + 




1 1121 

n 


1 + 


(1 + 


C21) 


Coi 1 
C02 ' 







(A19) 



and b) C02/C01 ^ 1 which then reduces the factor in the curly brackets to unity. These 
conditions ensure that the non-thermalized J=2 level has negligible interaction with the 
lower two either radiatively or coUisionally. Condition (a) is usually satisfied for optically 
thin mm/sub-mm lines but condition (b) imposes a temperature constraint, namely 
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In the case of the optically thin CI lines and since usually Tk ^ Tb ~ 2.7 K we use the 
neghgible radiation field approximation (f2i ~ fio ~ 1) to estimate the Qio and Q21 factors, 
which are then are given by 
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Using the detailed balance equations these become 
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In the weak radiation field limit the expression for K becomes 



(A21) 
(A22) 

(A23) 
(A24) 



K = 1 



1 



C20 
C21 



1121 _|_ Cio 
n C02 



1 



1 + ^ 
n 



+ 



C 



12 



1 + 



Cqi 

Cio 



1 + 



Coi 
nio 



l + Ci2/Cio \ nio 
,1 + Cio/CoJ n 
Coi 



1 + 



Cio 

C21 V C12/ V n / Coi V C20. 
or after some rearrangement of the terms and use of the detailed balance relations 



(A25) 



K = 1 + f 1 + ^) 1^ + ^e'^^/^^^ ( 1 + ^e-^^/^A ^G(n, Tk) f 1 + ^ 
V C21/ [ n g2 V go / C20 V n 

gl V §2 C21 J \ n 



nio \ ^ so^Ei/kTk I ^ ^ ai ^10^-Ei/kT, 
gl V go C20 



where we have set 
G(n,Tk) = l + ^e^^/'^^^ 



1 I 12 C2ig-E2i/kTk 
gl Cio 




(A26) 



(A27) 



-1 \ 1 + f^e^iATk 

The CI J=l-0, J=2-l hyperfine lines have degeneracy factors given by gj = 2J + 1, 
and energy levels Ei/k = 23.6 K, E2/k = 62.4 K. Their Einstein coefficients are Aio = 
7.93 X 10-^s-\ A21 = 2.68 x 10"V\ A20 = 2 x 10-^^s~\ and their collisional rates 710 = 



1.3 X 10"^° cm^ s-\ 721 



7.8 X 10-11 cm^ s-\ 720 



2.0 X 10 1° cm"* s 1 with a very 



weak temperature dependance (Zmuidzinas et al. 1988 and references therein). Thus the 
critical densities as defined previously are n2i = 964 cm~^ and nio = Aio/(7io + 712) = 
nio(O) [1 + (g2/gi)e"^2^/''i'''72i/7io]~^ The latter varies from ~ 600 cm^^ (Tk = K) to 
~ 300 cm~^ (high temperature limit). 

© 2004 RAS, MNRAS 000, [l]-?? 



CI lines as H2 gas mass tracers 27 

REFERENCES 

Aalto S., Booth R. S., Black J. M., & Johansson L. E. B. 1995 A&A, 300, 369 
Arimoto N., Sofue Y., & Tsujimoto T., 1996, PASJ, 48, 275 

Barone L. T., Heithausen A., Huttemeister S., Fritz T., & Klein U., 2000, MNRAS, 317, 649 

Baker A. J., Tacconi L. J., Genzel R., Lehnert M. D., & Lutz D. 2003 in The Neutral ISM in Galaxies, Marstrand, Sweeden, 
(in preparation) 

Barvainis R., Maloney P., Antonucci R., & AUoin D. 1997, ApJ, 484, 695 
Bash F. N., Green E., & Peters W. L. 1977, ApJ, 217, 464 
Bennett, G. L., et al. 1994, ApJ, 434, 587 
Blain A. W. 1999, MNRAS, 304, 669 

Blain A. W., Frayer D. T., Bock J. J., & Scoville N. Z. 2000, MNRAS, 313, 559 
Blitz L., & Shu F. H. 1980, ApJ, 238, 148 

Bolatto A. D., Jackson J. M., & Ingalls J. G. 1999, ApJ, 513, 275 

Bolatto A. D., Jackson J. M., Kraemer K. E., & Zhang X. 2000a, ApJ, 541, L17 

Bolatto A. D., Jackson J. M., Wilson G. D., & Moriarty-Schieven G. 2000b, ApJ, 532, 909 

Braine J., & Gombes F., 1992, A&A, 264, 433 

Bryant P. M., & Scoville N. Z., 1996, ApJ, 457, 678 

Brown R. L., & vanden Bout P. A. 1991, AJ, 102, 1956 

Brown R. L., & Vanden Bout P. A., 1992, ApJ, 397, Lll 

Garilli G. L., Menten K. M., & Yun M. S. 1999, ApJ, 521, L25 

Garilli G. L., Gox P., Omont A., Djorgovski S. G., Petric A., Beelen A., Isaak K. G., & McMahon R. G. 2002a, ApJ, 575, 145 

Garilli G. L. et al. 2002b, AJ, 123, 1838 

Garilli G. L., & Blain A. W. 2002c, ApJ, 569, 605 

Ghapman S. G., Blain A. W., Ivison R. J., & Small I. R. 2003, Nature, Vol. 422, pg. 695 
Ghiar J. E., Kutner M. L., Verier F., & Leous J., 1994, ApJ, 431, 658 
des Forets G. P., Roueff E., & Flower D. R. 1992, MNRAS, 258, 45 
de Boisanger G., & Ghieze J. P. 1991, A&A, 241, 581 

Devereux N., Taniguchi, Y., Sanders D. B., Nakai N., Young J. S., 1994, AJ, 107 (6), 2006 
Downes, D., Solomon, P. M. 1998, ApJ, 507, 615 

Downes D., Radford S. J. E., Greve A., Thum G., Solomon P. M., & Wink J. E. 1992, ApJ, 398, L25 

Dickman R. L., Snell R. L., & Schloerb F. P., 1986, ApJ, 309, 326 

Digel S., de Geus E., & Thaddeus P., 1994, 422, 92 

Ellis R. 2001, PASP, 113, 515 

Elmegreen B. G. 1989, ApJ, 338, 178 

Elmegreen B. G. 2002, ApJ, 577, 206 

Evans A. S., Frayer D. T., Surace J. A., & Sanders D. B. 2001, AJ, 121, 3285 

Falgarone E., Panis J.-F., Heithausen A., Stutzki J., Puget J.-L., & Bensch F. 1998, A&A, 331, 669 
Flower D. R., Le Bourlot J., des Forets G. P., & Roueff E. 1994, A&A, 282, 225 

Frayer D. T., & Ivison R. J., Scoville N. Z., Evans A. S., Yun M. S., Small I., Barger A. J., Blain A. W., & Kneib J.-P. 1999, 
ApJ, 514, L13 

Fixsen D. J., Bennett G. L., & Mather J. G. 1999, ApJ, 526, 207 

Genzel R., Baker A. J., Tacconi L. J., Lutz D., Gox P., Guilloteau S., & Omont A. 2003, ApJ, 584, 633 
Gerin M., & Phillips T. G., 2000, ApJ, 537, 644 

Greve T., Ivison R. J. I., & Papadopoulos P. P. 2003, ApJ (in press) 
Harrison A., Puxley P., Russel A., & Brand P. 1995, MNRAS, 277, 413 

© 2004 RAS, MNRAS 000, [l]-?? 



28 Papadopoulos, Thi, and Viti 

Hasegawa T., 1997, in lAU Symposium No. 170, CO: Twenty-five years of millimeter-wave spectroscopy, p. 39 (Kluwer Academic 
Publishers) 

Tielens A. G. G. M., & HoUenbach D. 1985a, ApJ, 291, 722 

Tielens A. G. G. M., & HoUenbach D. 1985b, ApJ, 291, 747 

HoUenbach D. J., Takahashi T., & Tielens A. G. G. M., 1991, ApJ, 377, 192 

HoUenbach D. J., & Tielens A. G. G. M. 1999, Reviews of Modern Physics, Vol 71, No 1, 174 

Ikeda M., Oka T., Tatematsu K., Sekimoto Y., Yamamoto S. 2002, ApJS, 139, 467 

Israel F. P., 1988, in Millimetre and Submillimetre Astronomy, R. D. Wolstencroft and W. B. Burton eds. Kluwer Academic 

Publishers, p. 281 
Israel F. P., 1997, A&A, 328, 471 

Israel F. P., Tilanus R. P. J., & Baas F. 1998, A&A, 339, 398 
Israel F. P., & Baas F. 2002, A&A, 383, 82 

Ivison R. J., Papadopoulos P. P., Seaquist E. R., & Bales S. A. 1996, MNRAS, 278, 669 
Jansen D. 1995, PhD Thesis, Leiden Observatory 

Johansson L. E. B., Olofsson H., Hjalmarson A., Gredel R., & Black J. H. 1994, A&A, 291, 89 
Jura M. 1975, ApJ, 197, 575 

Kaufman M. J., Wolfire M. G., HoUenbach D. J., & Luhman M. L., 1999, ApJ, 527, 799 
Keene J., Blake G. A., Phillips T. G., Huggins P. J., & Beichman C. A. 1985, ApJ, 299, 967 

Keene J., Lis D. C., Phillips T. G., & Schilke P. 1997, in Molecules in astrophysics: probes & processes lAU Symposium No 

178, ed. Ewine van Dishoeck (Cambridge: Cambridge Univ. Press) 129 
Langer W. D., & Penzias A. A. 1993, ApJ, 408, 539 

Lee H.-H., Herbst E., des Forets G. P., Roueff E., & Le Bourlot J. 1996, A&A, 311, 690 
Le Bourlot J., des Forets G. P., & Roueff E. 1993, ApJ, 416, L87 
Lioure A., & Chieze J. P. 1990, A&A, 235, 379 

MacLow M.-M., Klessen R. S., Burkert A., & Smith M .D. 1998, Phys. Rev. Lett., 80, 2754 

Madden S. C, PogUtsch A., Geis N., Stacey G. J., & Townes C. H. 1997, ApJ, 483, 200 

McCaU B. J., et al. 2003, Nature, 422, 500 

Maloney P. M., & Black J. H., 1988, ApJ, 325, 389 

Mehlert D., et al. 2002, A&A, 393, 809 

Meixner M., & Tielens A. G. G. M. 1993, ApJ, 405, 216 

Nakai N., & Kuno N. 1995, PASJ, 47, 761 

Neri R., Genzel R., Ivison R. J., Bertoldi F., Blain A. W., Chapman S. C, Cox P., Greve T. R., Omont A., & Frayer D. T. 
2003, |astro-ph/0307310| 

Ojha R., Stark A. A., Hsieh H. H., Lane A. P., Chamberlin R. A., Bania T. M., Bolatto A. D., Jackson J. M., & Wright G. A. 
2001, ApJ, 548, 253 

Oka T., Hasegawa T., Handa T., Hayashi M., & Sakamoto S. 1996, ApJ, 460, 334 
Oka T., et al., 2001, ApJ, 558, 176 

Omont A., Petitjean P., Guilloteau S., McMahon R. G., Solomon P. M., & Pecontal E. 1996, Nature, Vol 382, 428 
Pak S., Jaffe D. T., van Dishoeck E. F., Johansson L. E. B., & Booth R. S. 1998, ApJ, 498, 735 
Papadopoulos P. P., & Seaquist E. R., 1998, ApJ, 492, 521 
Papadopoulos P. P., & Ivison R. J. I., 2002, ApJ, 564, L9 
Papadopoulos P. P. 2004 (in preperation). 

Parravano A., HoUenbach D. J., McKee C. F. 2003, ApJ, 584, 797 

Pettini M., Steidel C. C, Adelberger K. L., Dickinson M., & Giavalisco M. 2000, ApJ, 528, 96 
Petitpas G. R., & Wilson C. D. 1998, ApJ, 503, 219 

© 2004 RAS, MNRAS 000, [l]-?? 



CI lines as H2 gas mass tracers 29 



Phillips T. G., & Huggins P. J. 1981, ApJ, 251, 533 
Rickard L. J., & Blitz L., 1985, ApJ, 292, L57 

Sakamoto S., Hasegawa T., Handa T., Hayashi M., & Oka T., 1997, ApJ, 486, 276 
Sanders D. B., Scoville N. Z., & Soifer B. T. 1991, ApJ, 370, 158 
Sawada T., et al. 2001, ApJS, 136, 189 

Schilke P., Carlstrom J. E., Keene J., & Phillips T. G. 1993, ApJ, 417, L67 

Solomon P. M., Scoville N. Z., & Sanders D. B., 1979, ApJ, 232, L89 

Solomon P. M., Downes D., & Radford S. J. E. 1992, ApJ, 398, L29 

Solomon P. M., Downes D., Radford S. J. E., & Barett J. W. 1997 ApJ, 478, 144 

Somerville R. S., Primack J. R., Faber S. M. 2001, MNRAS, 320, 504 

Spaans M., & van Dishoeck E. F., 1997, A&A, 323, 953 

Stark A. A., Bolatto A. D., Chamberlin R. A., Lane A. P., Bania T. M., Jackson J. M., & Lo K.-Y. 1997, ApJ, 480, L59 
Steidel C. C., Giavalisco M., Pettini M., Dickinson M., & Adelberger K. L. 1996, ApJ, 462, L17 
Stone J. M., Ostriker E. C., & Gammie C. F. 1998, ApJ, 508, L99 
Storzer, H., Stutzki J., & Sternberg A. 1997, A&A, 323, L13 
Stutzki J., et al. 1997, ApJ, 477, L33 

Suchkov A., Allen R. J., & Heckman T. M., 1993, ApJ, 413, 542 
Tauber J. A. 1996, A&A, 315, 591 

van Dishoeck et al. 1992, in lAU Symposium Astrochemistry of Cosmic Phenomena, P. D. Singh (ed.) pg 285 
Weiss A., Henkel C., Downes D., & Walter F. 2003, A&A, 409, L41 
Wilson C. D. 1997, ApJ, 487, L49 

Wilson C. D., Howe J. E., & Balogh M. L., 1999, ApJ, 517, 174 

White G. J., Ellison B., Claude S., Dent W. R. F., & Matheson D. N. 1994, A&A, 284, L23 
Wolfire M. G., McKee C. F., Hollenbach D., & Tielens A. G. G. M. 2003, ApJ, 587, 278 
Xie T., Allen M., & Langer W. D. 1995, ApJ, 440, 674 

Xie T. 1997, in CO: Twenty-five years of millimeter-wave spectroscopy, lAU Symposium No 170, W. B. Latter et al. (eds.) 

(Kluwer Academic Publishers) p. 131 
Yamamoto et al. 2001, ApJ, 547, L165 

Young J. S., & Scoville N. Z. S., 1991, ARAAS, Vol 29, 581 

Zhang X., Lee Y., Bolatto A., & Stark A. A. 2001, ApJ, 553, 274 

Zmuidzinas J., Betz A. L., Boreiko R. T., & Goldhaber D. M. 1998, ApJ, 335, 785 



© 2004 RAS, MNRAS 000, [T]-?? 



